Traveling waves in magnetized Taylor-Couette flow 
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We investigate numerically a traveling wave pattern observed in experimental magnetized Taylor- 
Couette flow at low magnetic Reynolds number. By accurately modeling viscous and magnetic 
boundaries in all directions, we reproduce the experimentally measured wave patterns and their 
amplitudes. Contrary to previous claims, the waves are shown to be transiently amplified distur- 
bances launched by viscous boundary layers rather than globally unstable magnetorotational modes. 
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I. INTRODUCTION 



The luminosity of most astrophysical accretion disks 
probably depends upon the magnetorotational instability 
(MRI) yy] , which has inspired searches for MRI in Taylor- 
Couette flow. Standard MRI modes will not grow un- 
less both the rotation period and the Alfven crossing 
time are shorter than the timescale for magnetic diffu- 
sion. This requires that both the magnetic Reynolds 
number Re^ = Q.iri{r2 — ri)lri and the Lundquist num- 
ber S = Vl^{r2 - ri)lri be > 1, where V^^ = B°/y/A^ 
is the Alfven speed and i?° is the background magnetic 
field parallel to the angular velocity fl = ile^. No labo- 
ratory study of standard MRI has been completed except 
for that of Sisan et al. [5], whose experiment proceeded 
from a background state that was already hydrodynami- 
cally turbulent before the field was applied. Recent linear 
analyses of axially periodic or infinite magnetized Taylor- 
Couette fiow has shown that MRI may grow at much re- 
duced i?ei„ and S in the presence of a combination of 
axial and current-free toroidal field [1, 0| ■ We call such 
modes "helical" MRI (HMRI) 

The PROMISE (Potsdam Rossendorf Magnetic Insta- 
bility Experiment) group have claimed to observe HMRI 
experimentally [1, [a, 01 ■ magnetic and flow parame- 
ters where linear analysis predicts instability, persistent 
fluctuations were measured that appeared to form axially 
traveling waves, consistent with expectations for HMRI. 
Similar behavior has been seen in nonlinear numerical 
simulations that approximate the experimental condi- 
tions, including realistic viscous boundary conditions for 
the velocities, but simplified ones for the magnetic field: 
perfectly conducting cylinders, and pseudo- vacuum con- 
ditions at the endcaps, when present 0, Both axially 
periodic and finite cylinders showed unsteady flow, the 
former case being more regular. However, the nonlin- 
ear simulations in [7|, l8| used somewhat different values 
for the cylinder rotation rates and other parameters than 
those reported in 



Previously, however, we have raised doubts about both 
the experimental realizability of HMRI and its astrophys- 
ical relevance [95. Finite cylinders with insulating endcaps 
were shown to reduce the growth rate and to stabilize 
highly resistive flows entirely, at least inviscid ones. 

Here we report nonlinear simulations with the ZEUS- 
MP 2.0 code [l^, which is a time-explicit, compress- 
ible, astrophysical ideal MHD parallel 3D code, to which 
we have added viscosity, resistivity (with subcycling to 
reduce the cost of the induction equation), and vac- 
uum boundary conditions, for axisymmetric flows in 
cylindrical coordinates (r, z)^]. The parameters of 
PROMISE as reported in or inferred from Q are used: 
gallium density p = 6.35 gcm~'^, magnetic diffusiv- 
ity 77 — 2.43 X 10^ cm^ s~^, magnetic Prandtl number 
= i^/rj — 1.40 X 10^^; Reynolds number Re = 
^iri{r2 — Tij/v = 1775; axial current 1^ = 6000 A; 
toroidal-coil currents I^p — 0, 50, 75, 120 A; and dimen- 
sions as in Fig. [1] For the flrst time, the flnite conduc- 
tivity and thickness of the copper vessel are allowed for 
('7Cu = 1.335 X 10^ cm^s~^), and this noticeably improves 
agreement with the measurements compared to previous 
linear calculations with radially perfectly conducting, ax- 
ially periodic boundaries Please note the difference 
of the direction of fJ, i?^ and (components measured 
in a right handed coordinate system) between this pa- 
per, where they are all assumed to be positive, and the 
experimental setup presented in [s^, where they are all 
negative (private communication). The direction of the 
traveling wave depends on the sign of the Poynting flux 
defined as —rriBipBz/ fJ^o Thus the direction of the 
traveling wave reported here is opposite as reported in 
5]. 



II. BOUNDARY CONDITIONS 



At the low frequencies relevant to PROMISE (/ 
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0.01 Hz), the skin depth of Copper 6^ = \/vCu/T^fl^o ~ 
19 cm, which is much larger than the thickness of the 
copper vessel surrounding the gallium in the PROMISE 
experiment, dw ~ 1.0 cm, so that the magnetic field dif- 
fuses rather easily into the boundary. On the other hand. 
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if one considers axial currents, the gallium and the copper 
wall act as resistors in parallel; taking into account their 
conductivities and radial thickness, one finds that their 
resistances are comparable [Rj : Ru : Rm =3:1:9; 
see Fig. [1] for the subscripts]. Therefore, the currents 
carried by the copper walls could be important for the 
toroidal field, and a perfectly insulating boundary condi- 
tion is also inappropriate. 

Figure 1: Computational domain for simulations of 
PROMISE experiment. Region (I): Inner copper cylinder, 
angular velocity Q\. (II): outer copper cylinder and bot- 
tom endcap, fl2- (HI): liquid gallium; (IV): vacuum. Thick 
dashed line: insulating upper endcap, = 0. Dimensions: 
ri = 4.0 cm; r2 = 8.0 cm; h = 40.0 cm; d^i ~ 1.0 cm; 
d^ii = 1.5cm; ni/2n = 3.6rpm; n2/2n = 0.972rpm. The 
exact configuration of the toroidal coils being unavailable to 
us, six coils (black rectangles) with dimensions as shown were 
used, with 67 turns in the two coils nearest the midplane and 
72 in the rest. Currents were adjusted to reproduce the 
reported fsl] Hartmann numbers Ha = B°ri / ^pjiQ-qu. 
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We have adapted a linear axisymmetric code developed 
by H, nil to allow for a helical field. Vertical periodicity is 
assumed, to allow separation of variables, but the full vis- 
cous and resistive radial equations are solved using finite 
differences, and a variety of radial boundary conditions 
can be imposed. For perfectly conducting boundaries and 
Itp = 75 A, where @| report persistent waves, our code in- 
deed finds a complex growth rate: s ~ 0. 0057-1-0. 057i s~^. 
But for insulating boundaries, the same parameters yield 
stability. 

This analysis points to the need for boundary condi- 
tions that accurately represent the influence of the copper 
vessels on the field. In the linear code just mentioned, we 
use the thin- wall approximation of [13l |. which in effect 
treats the cylinders as insulating for the poloidal field but 
conducting for the toroidal field. The errors of this ap- 



proximation increase with the ratio of wall thickness to 
gap width, which is not very small (« 0.25) in our case. 
Growth is predicted, but at a smaller rate than for per- 
fectly conducting walls, s « 0.0052 4- 0.056 z s""'^. The in- 
sensitivity of the imaginary part to the magnetic bound- 
aries supports the interpretation that these modes are 
hydrodynamic inertial oscillations weakly destabilized by 
the helical field [9]. 

In our nonlinear simulations, we include the copper 
walls (regions I and II) in the computational domain 
(Fig. [1]), but not the external coils themselves, whose in- 
ductive effects are therefore neglected. Outside the walls 
(region IV) we match onto a vacuum field Bcxt = V*!" 
vanishing at infinity. This is relatively straightforward 
in spherical geometry (used by many geodynamo exper- 
iments) because Laplace's equation separates. Our case 
is more difficult, because while Laplace's equation sep- 
arates in cylindrical coordinates when the boundary is 
an infinite cylinder, it does not fully separate outside a 
finite cylinder. Therefore we use an integral formula- 
tion that does not assume separability. The idea, called 
von Hagenow's method [13], is to find a surface current 
on the boundary that is equivalent to the current density 
in the interior as the source for Soxt via the free-space 
Green's function. The surface current is obtained by first 
solving the Grad-Shafranov equation [l^, [l^ in the in- 
terior with conducting boundary conditions, a problem 
that is separable in our case and is solved efficiently by 
combining FFTs along z with tridiagonal matrix inver- 
sion along r. 



III. RESULTS AND DISCUSSION 

We start with purely hydrodynamic (unmagnetized) 
simulations. For fx = 02/^^1 = 0.27, what we see is sim- 
ply an Ekman flow driven by the top and bottom end 
plates. Due to the stronger pumping at the upper, sta- 
tionary lid, the two Ekman cells are of unequal size. They 
are separated vertically by a narrow, intense radial out- 
flow, hereafter the "jet", lying at about 11cm above the 
bottom endcap. As discussed in [l3], the jet is unsteady 
at Re > 10^; it flaps or wanders rapidly in the poloidal 
plane. This has been verified by the PROMISE group 
(private communication). The amplitude of the flapping 
is ±0.4 mm s^i. 

Background states with purely axial or purely az- 
imuthal magnetic fields are symmetric under reflection 
z — i- — z, but a helical field breaks this symmetryfTs']. As 
a result, growing modes in vertically infinite or periodic 
cylinders propagate axially in a unique direction: that of 
the background Foynting flux —rflBgBz/ fio Fig. [2] 
displays vertical velocities near the outer cylinder in sim- 
ulations corresponding to the experimental runs of Q for 
several values of the toroidal current, Iip. A wave pattern 
very similar to that in the experimental data is seen. It 
is most obvious for = 75 A, just as in the experi- 
ment. Considering that we do not use exactly the same 
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Figure 2: (color). Axial velocities [mm s~^] versus time and 
depth sampled 1.5 cm Irom the outer cylinder, for the param- 
eters of the PROMISE experiment with toroidal currents 
as marked. Note height increases upward from the bottom 
endcap. No-slip velocity boundary conditions are imposed 
at the rigidly rotating endcaps, but the steady part of the 
resulting Ekman circulation is suppressed in these plots by 
subtracting the time average at each height. The waves ap- 
pear to be absorbed near the Ekman jet, at z ~ 100 mm. 
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Figure 3: (color). Panel (a): An extended version of the case 
— 75 A shown in Fig. [2] but without subtraction of the 
time average. The two Ekman cells are visible as the upflow 
(orange) at z > 100 mm and downflow (blue) at 2: < 100 mm; 
these are the expected directions of flow near the outer cylin- 
der. Panel (b): The same case again, except that after 
t = 360 s, the no-slip boundary condition at both endcaps is 
changed to an ideal Couette profile, i.e. 0,{r) = a + br~^ with 
a and b chosen to make 57 continuous at both cylinders; this 
eliminates Ekman circulation. Thereafter, the wave seems to 
be absorbed near the bottom (z ~ mm) rather than the jet 
{z ~ 100 mm), which itself dies out after t « 395 s. 
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external coil configuration as PROMISE, the agreement 
is remarkably good (Table. |T|. 



Table 1: Comparison of results for the frequency, wavelength, 
axial phase speed, and amplitude obtained from simulation 
and experiment for the case = 75 A. /i = f2i/27r is rotation 
frequency of inner cylinder. 
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^ 0.4 


^ 0.6 



Interestingly, the jet becomes nearly steady when I^p > 
50 A. It is known that Ekman circulation is significantly 
modified when the Elsasser number A = / {jiQprjVl) > 
If we use |B(ri)| for the field strength and ^2 for 
Q in this expression, then A = 4.8 at I,^ ~ 75 A. 

On the other hand, the magnetic field clearly promotes 
unsteadiness in the interior flow. The waves seen in Fig. [2] 
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are probably related to HMRI, but we do not believe that 
they arise from a global instability of the experimental 
Couette flow. To demonstrate this, we have repeated the 
third [I^p — 75 A) simulation shown in Fig. with differ- 
ent velocity boundary conditions. First, when we replace 
the rigidly rotating endcaps with differentially rotating 
ones that follow the ideal angular velocity profile of an 
infinitely long Taylor-Couette flow, then instead of the 
persistent travelling waves seen in Fig. [2l we see slowly 
damping standing waves, which we interpret as inertial 
oscillations excited by a small numerical force imbalance 
in the inital conditions Second, we perform a simula- 
tion that begins with the experimental boundary condi- 
tions until the traveling waves are well established, and 
then switches abruptly to ideal-Couette endcaps. After 
the switch, the Ekman circulation stops and the traveling 
waves disappear after one axial propagation time, as if 
they had been emitted by the Ekman layer at the upper 
endcap or by the layers on the upper part of the cylin- 



ders (Fig. [3]). After the switch in boundary conditions 
but before the waves fully disappear, their vertical phase 
speed increases from —0.7cm s~^ to — l.f cm s^^; the lat- 
ter is the speed predicted by linear analysis for axially 
periodic flow ^] (Fig. [3]). Both numerical tests support 
the interpretation that the wave pattern observed in the 
simulation and in the experiment is not a global HMRI 
mode but rather a transient disturbance that is somehow 
excited by the Ekman circulation and then transiently 
amplified as it propagates along the background axial 
Poynting flux, but is then absorbed once it reaches the 
jet or the bottom end cap. 
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